Transcriptome organization of white blood cells through gene co-expression network analysis in a large RNA-seq dataset.

Gene co-expression network analysis enables identification of biologically meaningful clusters of co-regulated genes (modules) in an unsupervised manner. We present here the largest study conducted thus far of co-expression networks in white blood cells (WBC) based on RNA-seq data from 624 individuals. We identify 41 modules, 13 of them related to specific immune-related functions and cell types (e.g. neutrophils, B and T cells, NK cells, and plasmacytoid dendritic cells); we highlight biologically relevant lncRNAs for each annotated module of co-expressed genes. We further characterize with unprecedented resolution the modules in T cell sub-types, through the availability of 95 immune phenotypes obtained by flow cytometry in the same individuals. This study provides novel insights into the transcriptional architecture of human leukocytes, showing how network analysis can advance our understanding of coding and non-coding gene interactions in immune system cells.

Hospital-associated venous thromboembolism prophylaxis use by risk assessment at a large integrated health care network in Northern California.

BACKGROUND: Hospital-acquired venous thromboembolism (HA VTE) is a preventable complication in hospitalized patients. OBJECTIVE: We aimed to examine the use of pharmacologic prophylaxis (pPPX) and compare two risk assessment methods for HA VTE: a retrospective electronic Padua Score (ePaduaKP) and admitting clinician's choice of risk within the admission orderset (low, moderate, or high). DESIGN, SETTINGS AND PARTICIPANTS: We retrospectively analyzed prophylaxis orders for adult medical admissions (2013-2019) at Kaiser Permanente Northern California, excluding surgical and ICU patients. INTERVENTION: ePaduaKP was calculated for all admissions. For a subset of these admissions, clinician-assigned HA VTE risk was extracted. MAIN OUTCOME AND MEASURES: Descriptive pPPX utilization rates between ePaduaKP and clinician-assigned risk as well as concordance between ePaduaKP and clinician-assigned risk. RESULTS: Among 849,059 encounters, 82.2% were classified as low risk by ePaduaKP, with 42.3% receiving pPPX. In the subset with clinician-assigned risk (608,512 encounters), low and high ePaduaKP encounters were classified as moderate risk in 87.5% and 92.0% of encounters, respectively. Overall, 56.7% of encounters with moderate clinician-assigned risk received pPPX, compared to 7.2% of encounters with low clinician-assigned risk. pPPX use occurred in a large portion of low ePaduaKP risk encounters. Clinicians frequently assigned moderate risk to encounters at admission irrespective of their ePaduaKP risk when retrospectively examined. We hypothesize that the current orderset design may have negatively influenced clinician-assigned risk choice as well as pPPX utilization. Future work should explore optimizing pPPX for high-risk patients only.

Echocardiographic heart ageing patterns predict cardiovascular and non-cardiovascular events and reflect biological age: the SardiNIA study.

AIMS: Age is a crucial risk factor for cardiovascular (CV) and non-CV diseases. As people age at different rates, the concept of biological age has been introduced as a personalized measure of functional deterioration. Associations of age with echocardiographic quantitative traits were analysed to assess different heart ageing rates and their ability to predict outcomes and reflect biological age. METHODS AND RESULTS: Associations of age with left ventricular mass, geometry, diastolic function, left atrial volume, and aortic root size were measured in 2614 healthy subjects. Based on the 95% two-sided tolerance intervals of each correlation, three discrete ageing trajectories were identified and categorized as 'slow', 'normal', and 'accelerated' heart ageing patterns. The primary endpoint included fatal and non-fatal CV events, and the secondary endpoint was a composite of CV and non-CV events and all-cause death. The phenotypic age of the heart (HeartPhAge) was estimated as a proxy of biological age. The slow ageing pattern was found in 8.7% of healthy participants, the normal pattern in 76.9%, and the accelerated pattern in 14.4%. Kaplan-Meier curves of the heart ageing patterns diverged significantly (P = 0.0001) for both primary and secondary endpoints, with the event rate being lowest in the slow, intermediate in the normal, and highest in the accelerated pattern. In the Cox proportional hazards model, heart ageing patterns predicted both primary (P = 0.01) and secondary (P = 0.03 to <0.0001) endpoints, independent of chronological age and risk factors. Compared with chronological age, HeartPhAge was 9 years younger in slow, 4 years older in accelerated (both P < 0.0001), and overlapping in normal ageing patterns. CONCLUSION: Standard Doppler echocardiography detects slow, normal, and accelerated heart ageing patterns. They predict CV and non-CV events, reflect biological age, and provide a new tool to calibrate prevention timing and intensity.

Age is the main risk factor for cardiovascular (CV) disease. Since people age and develop diseases at very different rates, biological age has been proposed as a more accurate measure of the body's functional decline. This study aimed to investigate the ageing rates of the heart and to assess their impact on CV events. The phenotypic age of the heart was also estimated as a proxy for biological age. Associations of age with Doppler echocardiographic parameters were analysed in a subgroup of 2614 clinically healthy subjects, part of a larger cohort of 3817 adults of both sexes.Three patterns of slow, normal, and accelerated ageing rates of the heart were detected. They predicted both CV and non-CV events, with different and progressively increasing event rates from the slow to the accelerated pattern. Compared with chronological age, the phenotypic (biological) age of the heart was 9 years younger in the slow pattern, 4 years older in the accelerated pattern, and comparable in the normal pattern.A standard Doppler echocardiogram is therefore able to detect three distinct heart ageing patterns, which reflect different biological susceptibilities to age-dependent diseases and provide a new tool for personalizing timeliness and intensity of prevention.

Glucose-6-phosphate dehydrogenase deficiency accelerates arterial aging in diabetes.

AIMS: High glucose levels and Glucose-6-Phosphate Dehydrogenase deficiency (G6PDd) have both tissue inflammatory effects. Here we determined whether G6PDd accelerates arterial aging (information linked stiffening) in diabetes. METHODS: Plasma glucose, interleukin 6 (IL6), and arterial stiffness (indexed as carotid-femoral Pulse Wave Velocity, PWV) and red blood cell G6PD activity were assessed in a large (4448) Sardinian population. RESULTS: Although high plasma glucose in diabetics, did not differ by G6DP status (178.2 ± 55.1 vs 169.0 ± 50.1 mg/dl) in G6DPd versus non-G6PDd subjects, respectively, IL6, and PWV (adjusted for age and glucose) were significantly increased in G6PDd vs non-G6PDd subjects (PWV, 8.0 ± 0.4 vs 7.2 ± 0.2 m/sec) and (IL6, 6.9 ± 5.0 vs 4.2 ± 3.0 pg/ml). In non-diabetics, neither fasting plasma glucose, nor IL6, nor PWV were impacted by G6PDd. CONCLUSION: G6PDd in diabetics is associated with increased inflammatory markers and accelerated arterial aging.

EDA ligand triggers plasma membrane trafficking of its receptor EDAR via PKA activation and SNAP23-containing complexes.

BACKGROUND: Ectodysplasin-A (EDA), a skin-specific TNF ligand, interacts with its membrane receptor EDAR to trigger EDA signaling in skin appendage formation. Gene mutations in EDA signaling cause Anhidrotic/Hypohidrotic Ectodermal Dysplasia (A/HED), which affects the formation of skin appendages including hair, teeth, and several exocrine glands. RESULTS: We report that EDA triggers the translocation of its receptor EDAR from a cytosolic compartment into the plasma membrane. We use protein affinity purification to show that upon EDA stimulation EDAR associates with SNAP23-STX6-VAMP1/2/3 vesicle trafficking complexes. We find that EDA-dependent PKA activation is critical for the association. Notably, either of two HED-linked EDAR mutations, T346M and R420W, prevents EDA-induced EDAR translocation; and both EDA-induced PKA activation and SNAP23 are required for Meibomian gland (MG) growth in a skin appendage model. CONCLUSIONS: Overall, in a novel regulatory mechanism, EDA increases plasma membrane translocation of its own receptor EDAR, augmenting EDA-EDAR signaling in skin appendage formation. Our findings also provide PKA and SNAP23 as potential targets for the intervention of HED.

Effect of Genetic Factors, Age and Sex on Levels of Circulating Extracellular Vesicles and Platelets.

Extracellular vesicles (EVs) mediate cell interactions in biological processes, such as receptor activation or molecule transfer. Estimates of variation by age and sex have been limited by small sample size, and no report has assessed the contribution of genetic factors to levels of EVs. Here, we evaluated blood levels of 25 EV and 3 platelet traits in 974 individuals (933 genotyped) and reported the first genome-wide association study (GWAS) on levels of these traits. EV levels all decreased with age, whereas the trend for their surface markers was more heterogeneous. Platelets and CD31dim platelet EVs significantly increased in females compared to males, although CD31 expression on both platelets and platelet EVs decreased in females. Levels of the other EV subsets were similar between sexes. GWAS revealed three statistically significant genetic signals associated with EV levels in the F10 and GBP1 genes and in the intergenic region between LRIG1 and KBTBD8. These add to a signal in the 3'UTR of RHOF associated with CD31 expression on platelets that was previously found to be associated with other platelet traits. These findings suggest that EV formation is not a simple, constant adjunct of metabolism but is under both age-related and genetic control that can be independent of the regulation of the levels of the cells from which the EVs derive.

Actively transcribed rDNA and distal junction (DJ) sequence are involved in association of NORs with nucleoli.

Although they are organelles without a limiting membrane, nucleoli have an exclusive structure, built upon the rDNA-rich acrocentric short arms of five human chromosomes (nucleolar organizer regions or NORs). This has raised the question: what are the structural features of a chromosome required for its inclusion in a nucleolus? Previous work has suggested that sequences adjacent to the tandemly repeated rDNA repeat units (DJ, distal junction sequence) may be involved, and we have extended such studies by addressing several issues related to the requirements for the association of NORs with nucleoli. We exploited both a set of somatic cell hybrids containing individual human acrocentric chromosomes and a set of Human Artificial Chromosomes (HACs) carrying different parts of a NOR, including an rDNA unit or DJ or PJ (proximal junction) sequence. Association of NORs with nucleoli was increased when constituent rDNA was transcribed and may be also affected by the status of heterochromatin blocks formed next to the rDNA arrays. Furthermore, our data suggest that a relatively small size DJ region, highly conserved in evolution, is also involved, along with the rDNA repeats, in the localization of p-arms of acrocentric chromosomes in nucleoli. Thus, we infer a cooperative action of rDNA sequence-stimulated by its activity-and sequences distal to rDNA contributing to incorporation into nucleoli. Analysis of NOR sequences also identified LncRNA_038958 in the DJ, a candidate transcript with the region of the suggested promoter that is located close to the DJ/rDNA boundary and contains CTCF binding sites. This LncRNA may affect RNA Polymerase I and/or nucleolar activity. Our findings provide the basis for future studies to determine which RNAs and proteins interact critically with NOR sequences to organize the higher-order structure of nucleoli and their function in normal cells and pathological states.

Development and Validation of a Diabetic Retinopathy Risk Stratification Algorithm.

OBJECTIVE: Although diabetic retinopathy is a leading cause of blindness worldwide, diabetes-related blindness can be prevented through effective screening, detection, and treatment of disease. The study goal was to develop risk stratification algorithms for the onset of retinal complications of diabetes, including proliferative diabetic retinopathy, referable retinopathy, and macular edema. RESEARCH DESIGN AND METHODS: Retrospective cohort analysis of patients from the Kaiser Permanente Northern California Diabetes Registry who had no evidence of diabetic retinopathy at a baseline diabetic retinopathy screening during 2008-2020 was performed. Machine learning and logistic regression prediction models for onset of proliferative diabetic retinopathy, diabetic macular edema, and referable retinopathy detected through routine screening were trained and internally validated. Model performance was assessed using area under the curve (AUC) metrics. RESULTS: The study cohort (N = 276,794) was 51.9% male and 42.1% White. Mean (±SD) age at baseline was 60.0 (±13.1) years. A machine learning XGBoost algorithm was effective in identifying patients who developed proliferative diabetic retinopathy (AUC 0.86; 95% CI, 0.86-0.87), diabetic macular edema (AUC 0.76; 95% CI, 0.75-0.77), and referable retinopathy (AUC 0.78; 95% CI, 0.78-0.79). Similar results were found using a simpler nine-covariate logistic regression model: proliferative diabetic retinopathy (AUC 0.82; 95% CI, 0.80-0.83), diabetic macular edema (AUC 0.73; 95% CI, 0.72-0.74), and referable retinopathy (AUC 0.75; 95% CI, 0.75-0.76). CONCLUSIONS: Relatively simple logistic regression models using nine readily available clinical variables can be used to rank order patients for onset of diabetic eye disease and thereby more efficiently prioritize and target screening for at risk patients.

A fast linkage method for population GWAS cohorts with related individuals.

Linkage analysis, a class of methods for detecting co-segregation of genomic segments and traits in families, was used to map disease-causing genes for decades before genotyping arrays and dense SNP genotyping enabled genome-wide association studies in population samples. Population samples often contain related individuals, but the segregation of alleles within families is rarely used because traditional linkage methods are computationally inefficient for larger datasets. Here, we describe Population Linkage, a novel application of Haseman-Elston regression as a method of moments estimator of variance components and their standard errors. We achieve additional computational efficiency by using modern methods for detection of IBD segments and variance component estimation, efficient preprocessing of input data, and minimizing redundant numerical calculations. We also refined variance component models to account for the biases in population-scale methods for IBD segment detection. We ran Population Linkage on four blood lipid traits in over 70,000 individuals from the HUNT and SardiNIA studies, successfully detecting 25 known genetic signals. One notable linkage signal that appeared in both was for low-density lipoprotein (LDL) cholesterol levels in the region near the gene APOE (LOD = 29.3, variance explained = 4.1%). This is the region where the missense variants rs7412 and rs429358, which together make up the ε2, ε3, and ε4 alleles each account for 2.4% and 0.8% of variation in circulating LDL cholesterol. Our results show the potential for linkage analysis and other large-scale applications of method of moments variance components estimation.

Polygenic risk score and biochemical/environmental variables predict a low-risk profile of age-related macular degeneration in Sardinia.

PURPOSE: To ascertain the prevalence and clinical and genetic features of age-related macular degeneration (AMD) in subjects living in the Lanusei valley, Central Sardinia, Italy, involved in a study on ageing (SardiNIA project). METHODS: A total of 814 volunteers aged ≥ 50 years, randomly selected from the SardiNIA project dataset, were included. A color fundus (CF) photograph of the 30° central retina of each eye was obtained and graded according to the Age-Related Eye Disease Study system. Life-style choices were investigated using standardized questionnaires. The concentrations of several inflammatory biomarkers (i.e., complement component, fibrinogen, and C-reactive protein) were measured. Polygenic risk score (PRS) was calculated and compared with results obtained from a European cohort. RESULTS: A total of 756 subjects had gradable CF photographs for AMD detection. In 91.3%, no signs of AMD were observed. The prevalence rates of early and late AMDs were 6.9% and 0.6%, respectively. A total of 85% of subjects were physically active; only 13.5% were current smokers. Low concentrations of complement component, fibrinogen, and C-reactive protein were found. We calculated the polygenic risk scores (PRS) using 40 AMD markers distributed on several candidate genes in Europeans and Sardinians. The mean PRS value was significantly lower in Sardinians than in the Europeans (0.21 vs. 0.248, respectively, p = 1.18 × 10-77). CONCLUSIONS: In our cohort, most subjects showed no sign of any AMD type and late AMD was a condition rarely observed. Results of genetic, biochemical, and life-style investigation support the hypothesis that Sardinia population may present of a peculiar background with a protective effect against AMD development.

Dilated hypertrophic phenotype of the carotid artery is associated with accelerated age-associated central arterial stiffening.

Hypertrophic carotid geometric phenotypes (h-CGP) are predictors of incident cardiovascular disease (CVD). While arterial aging is hypothesized as a contributor to this associated risk, the association of CGPs with chronological age is not clear. In this manuscript we examine whether hypertrophic CGPs represent accelerated biological, rather than chronological, aging by examining their association with carotid-femoral pulse wave velocity (PWV), the hallmark of arterial aging. We analyzed data from 5516 participants of the SardiNIA study with a wide range of age at baseline (20-101 years), and a median follow-up time of 13 years (mean 11.5 years; maximum 17.9 years). Baseline CGPs were defined based on the common carotid lumen diameter, wall thickness, and their ratio. Subject-specific rates of change of PWV, blood pressure parameters, body mass index, glucose, and lipids were estimated using linear mixed effects models. Compared to those with typical(t-) CGP, those with dilated hypertrophy (dh-) CGP had a greater longitudinal increase in PWV; this increase was significantly greater among older individuals and men. The greater PWV longitudinal increase in dh-CGP remained significant after adjusting for baseline values and rates of change of covariates. Dilated hypertrophic CGP is independently associated with accelerated increase in age-associated arterial stiffening over time, with a strong association in men than in women. Future studies are needed to examine if this association mediates the increased risk for CVD observed in individuals with hypertrophic cardiac remodelling and the role of retarding it to reduce this risk. HIGHLIGHTS: • Individuals with dilated hypertrophic geometric phenotypes of the common carotid artery (increased age- and sex-specific wall thickness and lumen diameter) have greater future central arterial stiffening, independently of other determinants of arterial stiffening. • The dilated hypertrophic phenotype group has a greater age-specific arterial dilation, wall thickening, and stiffness (the arterial aging triad). This suggests that this phenotype is a form of accelerated aging that might explain the worse clinic outcomes observed in this group. • Understanding the natural history of the carotid geometric phenotype across the lifespan and the determinants of the deleterious progression towards the dilated hypertrophic phenotype are needed to develop interventions that reduce the adverse clinical outcomes associated with it.

GWAS of genetic factors affecting white blood cell morphological parameters in Sardinians uncovers influence of chromosome 11 innate immunity gene cluster on eosinophil morphology.

Few genome-wide association studies (GWAS) analyzing genetic regulation of morphological traits of white blood cells have been reported. We carried out a GWAS of 12 morphological traits in 869 individuals from the general population of Sardinia, Italy. These traits, included measures of cell volume, conductivity and light scatter in four white-cell populations (eosinophils, lymphocytes, monocytes, neutrophils). This analysis yielded seven statistically significant signals, four of which were novel (four novel, PRG2, P2RX3, two of CDK6). Five signals were replicated in the independent INTERVAL cohort of 11 822 individuals. The most interesting signal with large effect size on eosinophil scatter (P-value = 8.33 x 10-32, beta = -1.651, se = 0.1351) falls within the innate immunity cluster on chromosome 11, and is located in the PRG2 gene. Computational analyses revealed that a rare, Sardinian-specific PRG2:p.Ser148Pro mutation modifies PRG2 amino acid contacts and protein dynamics in a manner that could possibly explain the changes observed in eosinophil morphology. Our discoveries shed light on genetics of morphological traits. For the first time, we describe such large effect size on eosinophils morphology that is relatively frequent in Sardinian population.

Genetic and biomarker modulation of arterial stiffness change in the SardiNIA population cohort.

Background and aims: Arterial stiffness (AS), quantified by pulse wave velocity (PWV), arises due to impaired arterial elastic tissue and smooth muscle dysfunction. We aimed to examine the longitudinal association of genetic, lipid and inflammation biomarkers with PWV and how these associations may change with aging. Materials and methods: We utilized genotype and four time-point biomarker data from the SardiNIA cohort [n = 6,301; mean baseline age 43.3 (SD 17.3); 58% females]. To investigate the association of PWV with genetic variants, lipid, and inflammation biomarkers, we employed linear mixed modeling, using age as the time scale. Biomarkers exhibiting significant longitudinal associations were categorized into tertiles and individuals within the second tertile or those with heterozygous alleles were excluded, leaving a cohort of 2,000 individuals. This cohort was further divided into four risk groups: low genetic and low biomarker (L-L), low genetic and high biomarker (L-H), high genetic and low biomarker (H-L), and high genetic and high biomarker risk (H-H). Subsequent analyses focused on these risk groups to assess their association to PWV with time. Results: Using the complete dataset, we found a significant longitudinal association of total cholesterol (TC), triglycerides (TG), fibrinogen (FGN), and total white blood cell count (TWBC) with PWV, all with p < 3.33 × 10-3. After grouping, individuals with homogeneous risk alleles of SNP rs3742207 and high baseline TG levels (H-H group) exhibited a 1.39-fold higher PWV (m/s) (95% CI, 1.17-1.64, p = 1.21 × 10-4) compared to the L-L group. Similarly, individuals in the H-H group of rs3742207-TWBC combination showed 1.75 times higher PWV (95% CI, 1.48-0.2.07, p = 1.01 × 10-10) compared to the L-L group. Similar patterns were observed for groups based on SNP rs7152623-TWBC risk. Furthermore, these associations became more pronounced with increasing age (p < 3.33 × 10-3). Conclusion: The longitudinal association of TG and TWBC biomarkers with PWV varied by SNPs rs3742207 and rs7152623 genotype. Further studies are warranted to investigate the function of genetics, lipids, and inflammation biomarkers on PWV change.

Personality traits are consistently associated with blood mitochondrial DNA copy number estimated from genome sequences in two genetic cohort studies.

Background: Mitochondrial DNA copy number (mtDNAcn) in tissues and blood can be altered in conditions like diabetes and major depression and may play a role in aging and longevity. However, little is known about the association between mtDNAcn and personality traits linked to emotional states, metabolic health, and longevity. This study tests the hypothesis that blood mtDNAcn is related to personality traits and mediates the association between personality and mortality. Methods: We assessed the big five personality domains and facets using the Revised NEO Personality Inventory (NEO-PI-R), assessed depressive symptoms with the Center for Epidemiologic Studies Depression Scale (CES-D), estimated mtDNAcn levels from whole-genome sequencing, and tracked mortality in participants from the Baltimore Longitudinal Study of Aging. Results were replicated in the SardiNIA Project. Results: We found that mtDNAcn was negatively associated with the Neuroticism domain and its facets and positively associated with facets from the other four domains. The direction and size of the effects were replicated in the SardiNIA cohort and were robust to adjustment for potential confounders in both samples. Consistent with the Neuroticism finding, higher depressive symptoms were associated with lower mtDNAcn. Finally, mtDNAcn mediated the association between personality and mortality risk. Conclusions: To our knowledge, this is the first study to show a replicable association between mtDNAcn and personality. Furthermore, the results support our hypothesis that mtDNAcn is a biomarker of the biological process that explains part of the association between personality and mortality. Funding: Support for this work was provided by the Intramural Research Program of the National Institute on Aging (Z01-AG000693, Z01-AG000970, and Z01-AG000949) and the National Institute of Neurological Disorders and Stroke of the National Institutes of Health. AT was also supported by the National Institute on Aging of the National Institutes of Health Grant R01AG068093.

Cells are powered by internal structures called mitochondria which have their own DNA molecules. How many copies of mitochondrial DNA blood cells contain is one aspect of mitochondrial health and is considered to provide a good indication of an individual's ability to convert glucose into energy. Consequently, changes in the amount of mitochondrial DNA in the blood are linked to conditions like diabetes and cancer, and have also been associated with aging and mortality. A set of well-classified personality traits known as 'the Big Five' have also been shown to affect energy levels and the longevity of individuals. However, it remained unclear if there is a relationship between these characteristics and the number of copies of mitochondrial DNA in the blood. To investigate, Oppong et al. used a specialized test to assess the personality traits of participants from two separate cohorts: Baltimore Longitudinal Study of Ageing and the SardiNIA Project. The genomic sequence of each person was then analyzed to calculate the amount of mitochondrial DNA in their blood, and their mortality was recorded based on whether they were alive or dead multiple years later. Oppong et al. found that low levels of mitochondrial DNA were linked with high scores in neuroticism (a trait typically associated with anxiety, depression, and self-doubt). Further statistical tests revealed that mitochondrial DNA levels mediate the relationship between a person's personality and their risk of death. These findings suggest that personality traits impact the number of mitochondrial DNA molecules in a person's blood, which, in turn, influences how long they are likely to live. However, further work is needed to find out what causes this effect.

Trends and Risk Factors for Venous Thromboembolism Among Hospitalized Medical Patients.

Importance: While hospital-associated venous thromboembolism (HA-VTE) is a known complication of hospitalization, contemporary incidence and outcomes data are scarce and methodologically contested. Objective: To define and validate an automated electronic health record (EHR)-based algorithm for retrospective detection of HA-VTE and examine contemporary HA-VTE incidence, previously reported risk factors, and outcomes. Design, Setting, and Participants: This cohort study was conducted using hospital admissions between January 1, 2013, and June 30, 2021, with follow-up until December 31, 2021. All medical (non-intensive care unit) admissions at an integrated health care delivery system with 21 hospitals in Northern California during the study period were included. Data were analyzed from January to June 2022. Exposures: Previously reported risk factors associated with HA-VTE and administration of pharmacological prophylaxis were evaluated as factors associated with HA-VTE. Main Outcomes and Measures: Yearly incidence rates and timing of HA-VTE events overall and by subtype (deep vein thrombosis, pulmonary embolism, both, or unknown), as well as readmissions and mortality rates. Results: Among 1 112 014 hospitalizations involving 529 492 patients (268 797 [50.8%] women; 75 238 Asian [14.2%], 52 697 Black [10.0%], 79 398 Hispanic [15.0%], and 307 439 non-Hispanic White [58.1%]; median [IQR] age, 67.0 [54.0-79.0] years), there were 13 843 HA-VTE events (1.2% of admissions) occurring in 10 410 patients (2.0%). HA-VTE events increased from 307 of 29 095 hospitalizations (1.1%) in the first quarter of 2013 to 551 of 33 729 hospitalizations (1.6%) in the first quarter of 2021. Among all HA-VTE events, 10 746 events (77.6%) were first noted after discharge. In multivariable analyses, several factors were associated with increased odds of HA-VTE, including active cancer (adjusted odds ratio [aOR], 1.96; 95% CI, 1.85-2.08), prior VTE (aOR, 1.71; 95% CI, 1.63-1.79), and reduced mobility (aOR, 1.63; 95% CI, 1.50-1.77). Factors associated with decreased likelihood of HA-VTE included Asian race (vs non-Hispanic White: aOR, 0.65; 95% CI, 0.61-0.69), current admission for suspected stroke (aOR, 0.73; 95% CI, 0.65-0.81), and Hispanic ethnicity (vs non-Hispanic White: aOR, 0.81; 95% CI, 0.77-0.86). HA-VTE events were associated with increased risk of readmission (hazard ratio [HR], 3.33; 95% CI, 3.25-3.41) and mortality (HR, 1.63; 95% CI, 1.57-1.70). Conclusions and Relevance: This study found that HA-VTE events occurred in 1.2% of medical admissions, increased over time, and were associated with increased adverse outcomes. These findings suggest that approaches designed to mitigate occurrence and outcomes associated with HA-VTE may remain needed.

Team is brain: leveraging EHR audit log data for new insights into acute care processes.

OBJECTIVE: To determine whether novel measures of contextual factors from multi-site electronic health record (EHR) audit log data can explain variation in clinical process outcomes. MATERIALS AND METHODS: We selected one widely-used process outcome: emergency department (ED)-based team time to deliver tissue plasminogen activator (tPA) to patients with acute ischemic stroke (AIS). We evaluated Epic audit log data (that tracks EHR user-interactions) for 3052 AIS patients aged 18+ who received tPA after presenting to an ED at three Northern California health systems (Stanford Health Care, UCSF Health, and Kaiser Permanente Northern California). Our primary outcome was door-to-needle time (DNT) and we assessed bivariate and multivariate relationships with six audit log-derived measures of treatment team busyness and prior team experience. RESULTS: Prior team experience was consistently associated with shorter DNT; teams with greater prior experience specifically on AIS cases had shorter DNT (minutes) across all sites: (Site 1: -94.73, 95% CI: -129.53 to 59.92; Site 2: -80.93, 95% CI: -130.43 to 31.43; Site 3: -42.95, 95% CI: -62.73 to 23.17). Teams with greater prior experience across all types of cases also had shorter DNT at two sites: (Site 1: -6.96, 95% CI: -14.56 to 0.65; Site 2: -19.16, 95% CI: -36.15 to 2.16; Site 3: -11.07, 95% CI: -17.39 to 4.74). Team busyness was not consistently associated with DNT across study sites. CONCLUSIONS: EHR audit log data offers a novel, scalable approach to measure key contextual factors relevant to clinical process outcomes across multiple sites. Audit log-based measures of team experience were associated with better process outcomes for AIS care, suggesting opportunities to study underlying mechanisms and improve care through deliberate training, team-building, and scheduling to maximize team experience.

Genetic Variation among Pharmacogenes in the Sardinian Population.

Pharmacogenetics (PGx) aims to identify the genetic factors that determine inter-individual differences in response to drug treatment maximizing efficacy while decreasing the risk of adverse events. Estimating the prevalence of PGx variants involved in drug response, is a critical preparatory step for large-scale implementation of a personalized medicine program in a target population. Here, we profiled pharmacogenetic variation in fourteen clinically relevant genes in a representative sample set of 1577 unrelated sequenced Sardinians, an ancient island population that accounts for genetic variation in Europe as a whole, and, at the same time is enriched in genetic variants that are very rare elsewhere. To this end, we used PGxPOP, a PGx allele caller based on the guidelines created by the Clinical Pharmacogenetics Implementation Consortium (CPIC), to identify the main phenotypes associated with the PGx alleles most represented in Sardinians. We estimated that 99.43% of Sardinian individuals might potentially respond atypically to at least one drug, that on average each individual is expected to have an abnormal response to about 17 drugs, and that for 27 drugs the fraction of the population at risk of atypical responses to therapy is more than 40%. Finally, we identified 174 pharmacogenetic variants for which the minor allele frequency was at least 10% higher among Sardinians as compared to other European populations, a fact that may contribute to substantial interpopulation variability in drug response phenotypes. This study provides baseline information for further large-scale pharmacogenomic investigations in the Sardinian population and underlines the importance of PGx characterization of diverse European populations, such as Sardinians.

Facets of Personality and Risk of Cognitive Impairment: Longitudinal Findings in a Rural Community from Sardinia.

BACKGROUND: Few studies have examined the associations between personality facets and dementia risk and rarely included individuals from rural settings or with low education. OBJECTIVE: To examine the association between personality and the risk of cognitive impairment. METHODS: Participants (N = 1,668; age 50 to 94 at baseline; 56.4% women; 86.5% less than high school diploma) were from a rural region of Sardinia (Italy) who completed the Revised NEO Personality Inventory (NEO-PI-R) during the first wave (2001-2004) and the Mini-Mental State Examination (MMSE) at waves two to five (2005-2021). Cox regression was used to test personality and covariates as predictors of cognitive impairment based on MMSE education-adjusted cutoffs. RESULTS: During the up to 18-year follow-up (M = 10.38; SD = 4.76), 187 individuals (11.2%) scored as cognitively impaired. Participants with higher neuroticism (particularly the depression facet [HR = 1.22, 95% CI = 1.06-1.40]), and lower agreeableness (particularly the modesty facet [HR = 0.83, 95% CI = 0.71-0.97]) and lower conscientiousness (particularly the dutifulness facet [HR = 0.78, 95% CI = 0.67-0.92]) were at higher risk of cognitive impairment. Lower warmth ([HR = 0.75, 95% CI = 0.65-0.87], facet of extraversion) and ideas ([HR = 0.76, 95% CI = 0.65-0.89], facet of openness) were also associated with increased risk of impairment. These associations were virtually unchanged in models that accounted for other risk factors, including smoking, depression, obesity, hypertension, diabetes, and apolipoprotein E (APOE) ɛ4 carrier status. Across the five domains, sex and the APOE variant did not moderate the associations. CONCLUSION: In a sample with demographic characteristics underrepresented in dementia research, this study identifies personality domains and facets most relevant to the risk of cognitive impairment.

Risk of severe COVID-19 infection among adults with prior exposure to children.

Susceptibility and severity of COVID-19 infection vary widely. Prior exposure to endemic coronaviruses, common in young children, may protect against SARS-CoV-2. We evaluated risk of severe COVID-19 among adults with and without exposure to young children in a large, integrated healthcare system. Adults with children 0-5 years were matched 1:1 to adults with children 6-11 years, 12-18 years, and those without children based upon a COVID-19 propensity score and risk factors for severe COVID-19. COVID-19 infections, hospitalizations, and need for intensive care unit (ICU) were assessed in 3,126,427 adults, of whom 24% (N = 743,814) had children 18 years or younger, and 8.8% (N = 274,316) had a youngest child 0-5 years. After 1:1 matching, propensity for COVID-19 infection and risk factors for severe COVID-19 were well balanced between groups. Rates of COVID-19 infection were slightly higher for adults with exposure to older children (incident risk ratio, 1.09, 95% confidence interval, [1.05-1.12] and IRR 1.09 [1.05-1.13] for adults with children 6-11 and 12-18, respectively), compared to those with children 0-5 years, although no difference in rates of COVID-19 illness requiring hospitalization or ICU admission was observed. However, adults without exposure to children had lower rates of COVID-19 infection (IRR 0.85, [0.83-0.87]) but significantly higher rates of COVID-19 hospitalization (IRR 1.49, [1.29-1.73]) and hospitalization requiring ICU admission (IRR 1.76, [1.19-2.58]) compared to those with children aged 0-5. In a large, real-world population, exposure to young children was associated with less severe COVID-19 illness. Endemic coronavirus cross-immunity may play a role in protection against severe COVID-19.

Genome-wide meta-analysis of iron status biomarkers and the effect of iron on all-cause mortality in HUNT.

Iron is essential for many biological processes, but iron levels must be tightly regulated to avoid harmful effects of both iron deficiency and overload. Here, we perform genome-wide association studies on four iron-related biomarkers (serum iron, serum ferritin, transferrin saturation, total iron-binding capacity) in the Trøndelag Health Study (HUNT), the Michigan Genomics Initiative (MGI), and the SardiNIA study, followed by their meta-analysis with publicly available summary statistics, analyzing up to 257,953 individuals. We identify 123 genetic loci associated with iron traits. Among 19 novel protein-altering variants, we observe a rare missense variant (rs367731784) in HUNT, which suggests a role for DNAJC13 in transferrin recycling. We further validate recently published results using genetic risk scores for each biomarker in HUNT (6% variance in serum iron explained) and present linear and non-linear Mendelian randomization analyses of the traits on all-cause mortality. We find evidence of a harmful effect of increased serum iron and transferrin saturation in linear analyses that estimate population-averaged effects. However, there was weak evidence of a protective effect of increasing serum iron at the very low end of its distribution. Our findings contribute to our understanding of the genes affecting iron status and its consequences on human health.